Strongly Localized State of a Photon at the Intersection of the Phase Slips in 2D 
Photonic Crystal with Low Contrast of Dielectric Constant 
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Two-dimensional photonic crystal with a rectangular symmetry and low contrast (< 1) of the 
dielectric constant is considered. We demonstrate that, despite the absence of a bandgap, strong 
localization of a photon can be achieved for certain "magic" geometries of a unit cell by introducing 
two 7r/2 phase slips along the major axes. Long- living photon mode is bound to the intersection 
of the phase slips. We calculate analytically the lifetime of this mode for the simplest geometry - 
a square lattice of cylinders of a radius, r. We find the magic radius, Tc, of a cylinder to be 43.10 
- - - percent of the lattice constant. For this value of r, the quality factor of the bound mode exceeds 

\ 10^. Small (~ 1%) deviation of r from Vc results in a drastic damping of the bound mode. 

' Introduction.-The field of photonic bandgap materials was pioneered by two papers In Ref. the idea of 

bandgap-induced inhibition of spontaneous emission, which limits the performance of light-emitting devices, was put 

forward. In Ref. |^ the observation was made that a photonic bandgap facilitates strong Anderson localization of 

light, which otherwise is hard to achieve jsj. Both papers addressed the issue of defect-induced localized photonic 

modes. However, the defects considered in Refs. and were of a different nature. In ||^ the defects were point-like 

scatterers, i.e. local perturbations of the spatially periodic dielectric constant. In contrast, a defect discussed in Q] 

represented a periodicity interruption or, in other words, a phase slip While both types of defects cause the 

I ■ formation of the in-gap states, there is a fundamental difference between them. A point-like defect does not lift the 

^ long-range order of the underlying crystal. On the contrary, in the presence of a phase-slip, the distance between two 

lattice sites located to the left and to the right from the phase slip is non- integer (in the units of the lattice constant). 

Lately, the properties of defects in photonic bandgap structures (phase-slip-like and more complex defects, 

specifically designed for certain manipulations of the light flow [Q) have become a branch of research of their own 

Theoretical studies of photonic bandgap materials are almost exclusively focused on structures with a high contrast 

of the dielectric constant which is sufficient for the formation of either 2D or 3D bandgap, i.e. the region of frequencies 

' where the propagation of light is completely forbidden. For low contrast, periodicity manifests itself in the formation 

"j--; of narrow stop bands only along certain directions close to the Bragg condition. The overall density of photonic states 

is weakly perturbed by such a periodicity. This rules out the possibility of defect-induced localized modes. Although, 

^ ! strictly speaking, this statement is correct, we demonstrate in the present paper that strongly localized photon modes 

I ' are possible in photonic crystals with contrast of dielectric constant < 1. By "strongly localized" modes we mean 

' O ' the modes with decay time much longer than the inverse frequency, or, in other words, with high quality factor, Q. 

^ ' In fact, we will demonstrate below that, when the contrast of the dielectric constant is low, the proper choice of the 

geometry and parameters of the 2D unit cell, combined with proper periodicity interruptions, allows to achieve a 

Q- factor of the locaHzed mode as high as Q ~ 10^. 

Qualitative Consideration - In the absence of photonic gap, it can be concluded on general grounds that formation 

of a defect-induced localized mode is impossible, since its frequency would be degenerate with the continuum of the 

^ ' extended Bloch waves. There is, however, an exception. Consider a 2D photonic crystal in which dielectric constant 

5^ ' e(x, y) is separable, e{x, y) = eo + ei(a^) + ^2{y), where ei, and £2 are weak periodic modulations with periods a and b, 

respectively. Thus, for the light wave propagating strictly along the a;-direction the position of a narrow stop band 

1/2 

in the dispersion law is determined by the Bragg condition Eq uJx — nc/a. Analogously, for the light wave along the 

1 /2 

y-direction, the Bragg condition reads Eq ujy — nc/b. Assume now, that the phase slips along both directions, x and 

y, are introduced into the photonic crystal. This situation is illustrated in Fig. 1. Note, that, even in the presence of 

the phase slips, e{x,y) remains separable. On the other hand, it is known that, much like a point defect, a phase slip 

in one dimension causes a localized photonic state that decays exponentially away from the phase slip Then, in 

the presence of two phase slips, a truly localized photonic mode bound to the intersection of the phase slips, x — y — 0, 

1/2 

exists. The frequency of this mode is close to luq — (w^ + ujy) . We would like to emphasize that existence of a 
2D localized mode, degenerate with continuum, is unique for the case of the phase slips, since this mode is effectively 
decoupled from the continuum. On the other hand, even for separable e{x,y), a point defect of arbitrary strength is 
unable to localize a photon when the gap is absent. This is because, in the presence of a defect, the separability is 
lifted. 

Naturally, any violation of separability of e{x, y) would result in the leakage of a mode bound to the intersection 
of the phase slips. It is also obvoius that separability is lacking for any realistic geometry of photonic crytal (like the 
one illustrated in Fig. 1). Our prime observation is that, for certain "magic" geometries of the unit cell, the leakage 
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can be suppressed even without separabihty. Then the Hfetime of the mode remains much larger than uJq^- 

We start with the remark that, unUke the case of a random disorder, the decay of the locahzed mode is, to the 
first approximation, determined by a single Fourier component of e{x,y). Our main idea is that this component 
can he eliminated by a proper arrangement of a unit cell. As the simplest example, consider a rectangular lattice of 
cylinders of radius r (Fig. la). In this case, the Fourier component, responsible for the leakage, is proportional to 

Ji ^27r-\/(r/a)2 + {r/by^ , where Ji is the Bcsscl function of the first order. Correspondingly, the magic value of r, 
for which the leakage is suppressed, is determined by a first zero of Ji{u), i.e. u — uq ^ 3.830. This yields for the 
magic radius, Tc, the value Tc = 0.610a6/\/a^ + In particular, for a square lattice, a 
Away from the magic radius, the quality factor falls off as (r — Tc)"^, 



6, we obtain Tr — 0.4310a. 
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where 5t is the difference between dielectric constants of the cylinder and the background. For the coefhcient C in 
the most interesting case of 7r/2 phase slips we obtain below 



C = 
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Ji(2-i/2«o) 
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where Jq is the Bessel function of a zero order. In the limit r ^ Tc the growth of Q saturates due to the higher-order 
processes illustrated in Fig. 2. These processes involve, strictly speaking, all the Fourier components with non-zero 

n and m; these components are proportional to Ji ^27r^(nr/a)^ -f (mr/b)^^ . Our important observation is that 

the leakage can be further suppressed by fine tuning of r around re- This is due to the compensation of different 
contributions - phenomenon analogous to the Fano resonance The final expression for the quality factor, which 
incorporates the higher-order processes, has the form 



-i = c 
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where dimensionless factors ki and K2 involve the following lattice sums 
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The elements Fn.m related to the Fourier components of e(x,?/) as follows 

The sum in Eq. (|^) reflects the processes 2 — 2i in Fig. 2, while the sum in Eq. (||) originates from the process 2 — 22. 
First term in Eq. (|) describes the Fano resonance. It suggests that the leakage can be reduced, if r slightly exceeds 
Tc by (Se/eo)Kia. Then the maximal possible value of the quality factor is given by 



Cn,2 



Formula (0) is our main quantitative result. In general, it could be expected that parameters ki and K2 are ^ 1. 
Remarkably, their values turn out to be very small, ki « 5 • 10~^ and K2 ~ 0.8 • 10^"^. Thus, even for Se — eo, we 
get a very high quality factor, Qm, ~ 0.4 • 10^. Achieving this value, however, requires a high precision in the choice 
of r. For example, without fine tuning (i.e. for r — rc), Eq. ^ suggets that the quality factor drops from Q,n to 

« 0.03Q™. 
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Derivation of Eqs. ^).- For simplicity we consider the case of TM polarization with electric field, £{x, y), along 
the z — axis. Since the derivation is based on the perturbation theory, it is convenient to cast the wave equation for 
£{x,y) into Schrodinger-like form 

where A = stands for the "energy", while the "potentials" Ui{x), U2{y), and U^Jl{x,y) are defined as 

Ui{x) = ^ Unfi{x + d;,;Sign(a;)), U2{y) = ^ UQ^rn{y + dySign{y)), (9) 

n m 



Upert{x,y)^ ^ U^P^{x,y)= Un,mix + d^sign{x),y + dySigii{y)). 



(10) 



n,m>0 



n,m>0 



The arguments in r.h.s. of Eqs. (^), ( p^ reflect the periodicity interruptions by 2d,j. and 2c?y, shown in Fig. 1(a), 
while the functions Un^m{x,y) are the components of the Fourier expansion of periodic e(x,y) in the absence of the 
phase slips 



(11) 



where gx = 7r/a and gy = n/b. In these notations, the known solution for the localized mode induced by a 
one-dimensional phase slip takes the form 



E^„''\x) - (27,)i/2cos(<?,x)exp(-7,|a:|), E^^\y) = {2jyy/^ cos{gyy) exp {-jy\y\) , 
where the decrements "fx, ly are given by 
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Therefore, the zero-order solution of Eq. (^) corresponding to U^^fl = is £Q{x,y) ~ e\^^ {x)e'^^ (y) with the "energy" 



Ao = gl + gl + 5e 



Fifl\ cos 



/ 2'Kdx 



Fo 1 1 cos 
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The contour of equal intensity for £(x,y) is shown schematically in Fig. lb. 

Calculation of leakage of the localized mode S^^x^y) requires taking into account the phasc-slip-induced perturba- 
tion, Up^^l{x,y), up to the second order. We will restrict further consideration to the most interesting case of 7r/2 
phase slip, i.e. dx = a/4, dy = b/A. For this purpose we will need, in addition to £a{x^ y), the zero-order solutions of 

Eq. (||), corresponding to the continuous spectrum of propagating modes. These solutions £^'^"(0;, y) = Elfl{x)Ei^q{y) 
can be easily found in a similar way as £o{x, y) 



V(i+^ 



|e '^"^1 sin(pa;)| - e*^"^ pcos{px) + ifj.{l + p^)^^^ sm{px) |, 



(15) 



where p = 2pgxeo/{Se \Fi^q\), and the index /z ± 1 enumerates the upper and the lower branches of a one-dimensional 
spectrum in the a;-direction; these branches are separated by a narrow Bragg gap Se \Fifl\ (see Fig. 2). The "energies" 
corresponding to £!^'^{x,y) are given by 



-p,9 - gl + gl + ^^ 



(16) 



The form of the wave equation (g) allows to write down automatically the epression for the first and second-order 
corrections, Aq^'', Aq^'', to the "energy" Aq, Eq. (|lj). 



3 



AO — Ap,q 



fi,p,l'.q 



Obviously, Aq^'' does not cause any leakage, and, moreover, Aq^'' = due to symmetry of the geometry under consid- 
eration. First non-vanishing contribution to ImA, which describes the leakage, as follows from the golden rule, 

ImA = ImA^') = ^ ^ K^o 1^1,2 |£;;,'')P<5(Ao - A^^-;), (18) 

emerges if the lowest order expression for the operator Ti,2 

f['l^u[f{x,y) + uif{x,y) (19) 

is substituted into Eq. (|l8|). Evaluation of ImA leads to Eq. (|l|) for the quality factor, Q = 2A/ImA = lo/Ivcluj. Since 
for magic configuration this term is zero, the operator T12 should be taken with higher accuracy. Namely, higher 
order processes, illustrated in Fig. 2, amount to the following modification of the operator Ti 2 

where Ep^^py = ex]){ipxX + ipyu), Xp^^p^ — Px + Py stand for "nonresonant" solutions of Eq. (^) with \p\ > l-nja. For 

these solutions, periodic modulation can be neglected. Matrix elements of C/n^m between E^ix^y) and plane waves, 
required for evaluation of Eq. ( [T^ ) , can be calculated analytically. The form of these matrix elements suggests that 
the main contribution to ImA comes from the small angular regions ~ (5e/eo in Fig- 2 (of the order of the Bragg gap) 
around the Bragg directions. After that, the summation in Eq. (|l8|) with the use of the 5-function, reduces to the 
integrals 

In= duu{l+ur^'-+''^/'= (21) 

Jo (2n+l)(2n + 3) 

with n = 1, 2, and 3. Now, as summation over all momenta in Eqs. ([l^), ( pO| ) is performed analytically, we arrive at 
Eq. (D. 

Conclusion- First photonic-crystal based lasers were reported in Refs. ]ll| , p^ . In Ref. a two-dimensional 
photonic crystal was fabricated within the active region of the InGaAsP multiple quantum well structure using the 
etching procedure. The contrast of the dielectric constants between the semiconductor and etched air-holes was high 
(« 10 : 1). This contrast was sufficient for opening of the two-dimensional bandgap. Then a defect, which was 
essentially a missing air-hole |Q , gave rise to the in-gap mode localized within approximately two periods of the host 
hexagonal crystal. It was demonstrated in Ref. ||ll[] that above the threshold the laser emission is dominated by this 
defect mode. 

A low-contrast photonic crystal with array of holes etched in the 5*102 substrate (e = 1.46) of a solid organic gain 
film (e = 1.7) was reported in Ref. |1^. Without photonic gap, the authors did not attempt to create a localized 
mode, although in the later work |1^ a square symmetry of the lattice pattern was employed. 

In the present paper we have demonstrated that a localized photon mode is possible in a 2D photonic crystal with a 
low contrast of the dielectric constant. In a sense, the absence of photonic gap can be compensated by a high-precision 
tailoring of the unit cell parameters. There is another example, when an accurate choice of point-defect parameters 
gives rise to a long lifetime. It pertains to the crystals with a high contrast of the dielectric constant jljjl^. Namely, 
in a finite thickness 2D photonic crystal with a complete gap for in-plane light propagation, the leakage of the in-gap 
mode in the ^-direction can be strongly suppressed for a certain radius of a cylinder constituting a defect . 

A natural extension of the present study is to explore a possibility to guide light from one location to another using 
an array of phase slips in low-contrast magic crystals. It would be appealing if such crystals offered an alternative to 
the conventional bandgap materials Q in terms of manipulating the light flow. 

This work was supported by the Army Research Office under Grant No. DAAD 19-0010406 and the Petroleum 
Research Fund under Grant No. 37890- AC6. 
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FIG. 1. (a) Schematic illustration of two intersecting phase slips in a 2D photonic crystal with rectangular symmetry, (b) 
The contour of equal intensity for the mode localized at the intersection. 
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FIG. 2. k-plane for a photonic crystal with low contrast of dielectric constant. Solid line is a contour of constant \k\ = Cq^^w/c. 
Stop bands near the Bragg conditions are indicated with long-dashed lines. Thick solid lines: (1) - second-order process 
responsible for the decay of the bound mode; (2) - fourth-order process responsible for the decay of the bound mode. Virtual 
transition 2 is followed by one of the virtual transitions 2i or 22 shown with dashed lines. 
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